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ESTIMATION OF DOPPLSR SHIFT COMPENSATION 
IN A MOBILE CQMMmEcATIONrTYSTEM 



The present invention relates to estimating a Doppler shift 
Gompensation in a mobile communication system. 

In a mobile communication system, signals which are transmitted 
from mobile stations moving relative to a base station are 
subject to the well-known , Doppler effect, which causes a 
frequency shift in the frequency received at the base station 
relative to that which was transmitted by the mobile station. 
This frequency shift is referred to herein as the Doppler shif t • 
The Doppler shift is dependent upon the speed and direction of 
movement of the mobile station relative to the base station. 
Thus, the Doppler effect can provide an increase or a decrease 
in the frequency, depending on the direction of movement of the 
mobile station relative to the base station. The magnitude of 
the Doppler shift is dependent on the speed with which the mobile 
station is moving relative to the base station. 

Another area of frequency offsets arises when the transmitter 
unit and the receiver unit are wrongly synchronised. 

Existing mobile communication installations provide a form of 
Doppler compensation, in that the frequency detection circuitry 
within the base station which selects a particular signal on a 
particular channel can take into account a certain amount of 
Doppler shift in the signal. 

AU 664626 relates to a method and circuit arrangement for 
compensating for the Doppler shift in a radio signal propagating 
between a base station and a mobile station when the mobile 
station approaches and moves past the base station. As a mobile 
station approaches the base station with a decreasing propagating 
time delay, the propagating time delay is integrated at intervals 
to determine how the propagating time delay varies with time. 
This information is used to effect a change in the radio signal 
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frequency at a particular time to compensate for the sudden 
Doppler shift as the mobile station moves past the base station. 
In a TDMA system, the interrogation occurs at intervals equal to 
an integral number of t ime frames* In the GSM standard a time 
frame comprises eight consecutive time slots and a single 
transmission burst passes between a particular mobile station and 
base station in any one time frame. 

Thus, in this method of Doppler compensation decisions are made 
in response to past and incoming signals to improve the reception 
at a future time; i.e. it is a reactive system. It would be 
desirable to implement a system which can actively compensate for 
Doppler shifts in an incoming signal in real time. 

Moreover, due to fading a signal subject to the Doppler effect 
may have a variation in amplitude as well as frequency. It would 
be desirable to estimate the Doppler effect on amplitude as well 
as frecpiency. 

According to one aspect of the invention there is provided a 
system for generating a Doppler correction factor for 
compensating for Doppler shift in a signal transmitted between 
a mobile station and a base station in a mobile communication 
system, the system comprising: a channel impulse response 
determination circuit for deteirmining a channel impulse response 
for the. channel on which the signal is received; reference 
circuitry for generating a reference vector from a set of known 
symbols and the channel impulse response; means for generating 
an error vector from the reference vector and received samples 
of the signal; means for determining noise variance of the 
incoming signal; and an estimator for generating a Doppler 
correction factor having a real part representing amplitude 
correction and an imaginary part representing phase correction, 
the estimator utilising the noise variance^ the error vector, the 
samples of the received signal and the sampling times. 

According to the GSM standard, a TDMA transmission burst 
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comprises a training sequence which is normally used to determine 
the channel impulse response for the channel on which the signal 
is received. This .is done by a convolution of the received 
training sequence with a stored version of the training sequence. 
Thus, the training sequence can constitute the set of known 
symbols for use in generating the reference vector. 
Alternatively, the set of known symbols can comprise estimated 
symbols from the received signal samples. 

In a TDMA system, the signal comprises a sequence of transmission 
bursts in respective time slots. The Doppler correction factor 
can be generated for each transmission burst and used to correct 
signal samples in that burst , 

In addition for correcting for the Doppler effect, the Doppler 
correction factor also corrects for frequency offsets arising 
from wrong frequency synchronisation between transmitter and 
receiver units- 

According to a further aspect of the present invention there is 
provided a method for compensating said Doppler shift in a signal 
transmitted between a mobile station and a base station in a 
mobile communication system, the method comprising: determining 
a channel impulse response for the channel on which the signal 
is received; generating a reference vector from a set of known 
symbols in the channel impulse response; generating an error 
vector from the reference vector and received samples of the 
signals; determining noise variance of the incoming signals; 
estimating a Doppler correction factor having a real part 
representing amplitude correction and an imaginary part 
representing phase correction, the estimating step using the 
noise variance, the error vector, the samples of the received 
signals and the sampling times; and applying the Doppler 
correction factor to the received signal samples thereby to 
compensate for Doppler shift. 

For a better understanding of the present invention and to show 
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how the same may be carried into effect, reference will now be 
made by way of example to the accompanying drawings in which: 

Figure .^^,,^1^ is a diagram of a signal burst in a mobile 
communication system; 

Figure 2 is a block diagram of circuitry for implementing 
Doppler shift compensation; and 

Figure ^ 3 _i s a diagram illustrating the Estim /3 block. 

Figure 1 illustrates a normal burst in a mobile communication 
system according to the GSM standard* This figure represents a 
burst received at a base station. For a TDIVIA system according 
to the GSM standard, mobile stations transmit bursts as modulated 
signals on frequency channels allocated by a base station 
controller. One frequency channel may support up to eight 
bursts, each burst associated with a respective call, where each 
call is allocated a time slot in which to send the burst. 
Further details of a TDMA system according to the GSM standard 
are not described herein because they are known to a person 
skilled in the art. 

The normal burst contains two packets of 58 bits (DATA) 
surrounding a training sequence (TRS) of 26 bits. Three tail 
bits (TS) are added at each end of the normal burst. The 
training sequence (TRS) is a predetermined sequence of bits which 
is sent by the mobile station (MS) and is known at the base 
station controller (BSC) . In Figure 1, it is assumed that a- 
received signal sample y^ represents the kt^ bit of the burst, in 
this instance the start of the training sequence. The training 
sequence is utilised at the base station controller to estimate 
the impulse response of the channel over which the burst is sent, 
and, at least in one embodiment of the invention to estimate a 
Doppler compensation factor. The actual information which is 
transmitted is located in the data bits (DATA) of the burst , 

As explained earlier, the environment through which a signal 
passes from a mobile station to a base station can vary 
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considerably, depending, amongst other things, on the distance 
between the mobile station and the base station, and interference 
caused by buildings and other structures in the area. As a 
result, the signal strength and signal quality of the signal 
received at the base station varies widely. Moreover, for moving 
mobile stations, the signal received by the base station is 
subject to a Doppler shift which should be corrected. 

Figure 2 illustrates a circuit 1 suitable for implementing a 
Doppler compensation in a GSM system. It should be understood 
that the various blocks in Figure 2, although illustrated as 
separate interconnected entities, do not necessarily represent 
separate physical entities, but are intended to represent 
diagrammatically the various steps which are carried out . The 
blocks could be implemented as circuits or a suitably programmed 
microprocessor may effect each of the functions which is 
individually assigned to the blocks. 

An antenna 20 receives signals 11 from the mobile stations MS. 
The antenna 20 is connected to RF circuitry 22. This circuitry 
22 operates on the received- burst to downshift the frequency to 
the baseband frequency and to sample the burst to provide from 
the analogue signal digital sampled values y^ where i is the 
sample index. The RF circuitry 22 outputs the sampled values 
yo- . .y^ for each burst together with the sampling times to^.-t^i, 
typically at the e2q>ected bit rate of the transmitted signal* 
The output of RF circuitry 22 is supplied along line 24 to a 
channel impulse response (C.I.R.) block 10, to a variance and 
error calculator 16 to filtering and equalisation circuitry 12, 
to an Estim {3 block 3 6 for estimating a Doppler correction factor 
and to transforming circuitry 4 0 to enable the application of the 
Doppler correction factor to the burst. 

A memory 32 holds the training sequence x which is the 
predetermined sequence of bits which is sent by the mobile 
station MS as a training sequence and received at the base 
station and converted into samples y^^ = Yi^-^y^+ss- "^^ae reference 
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training sequence x is supplied to a reference generator 14 and 
to the channel impulse response (C.I.R.) block 10, The reference 
generator 14 also receives the estimated channel impulse response 
h from the channel impulse response block 10 . 

The C.I.R. block 10 receives the burst, including the received 
training sequence y^^^ ^'^^ calculates an estimated channel impulse 
response h for example by calculating the cross correlation 
between the received training sequence y^Rs ^^<^ t.he known training 
sequence x according to equation 1. There are other types of 
estimation methods which could be used. 

h = xcorr (x^y^^^) (equation 1) 

It will be appreciated that, prior to effecting the cross 
correlation, the known training sequence x, which is stored in 
digital form, is i,q modulated similarly to the manner in which 
the training sequence has been modulated at the MS for 
transmission, according to the GSM standard. The cross 
correlation is done in a known manner to produce a channel 
impulse response in the form of five tap values (h(j)j.oto4) • 

As is known, the estimated impulse response h is used to 
calculate the expected estimate of the data in the received 
burst, as though the data has been subject to the same average 
noise . 

The CI.R. block also generates timing advance information t, 
which is used to determine where in the allocated time slot the 
received burst is located. 

For each burst, the estimated channel impulse response h for that 
burst is calculated by the CIR block 10 and is supplied to 
filtering/equalising circuitry which allows the data, DATA(r) , 
in that burst to be recovered. As is known, the 

filtering/equalising circuit 12 receives the channel impulse 
response h and timing information r for the received burst to 
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allow the signal to be demodulated, filtered and decoded, to 
recover the data in a known manner. 

The reference generator 14 produces a reference vector, z, which 
is calculated using the convolution of the impulse response and 
the known training sequence . The reference generator 14 performs 
the following calculation: 

z = h*x (equation 2) 

where x represents the known symbol values depending o the used 
modulation . 

z„ represents the kth sample of the vector z and z denotes a 
diagonal matrix of the vector samples Zj,, that is, 2=xh, where x 
is the diagonal matrix of the known training sequence symbols. 

The reference vector z is supplied from the reference generator 
to the variance and error calculator 16. As described above, the 
variance calculator also receives the samples Yk including the 
received training sequence y^s • The variance calculator 
calculates a variance var (a^) according to the following 
equation : 

25 

var = EdYi - Zi I '^2)) (equation 4) 

i=4 



ref_length 

The term ref_length is a constant representing the length of the 
reference vector, z. This is calculated by multiplying the 
number of samples (22) by the bit separation. 

In equation 4, the values of Yl are the sampled values of the 
received training sequence for the burst (yk-*yk+25) • 



The variance and error calculator also generates an error vector 
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z-y {i.e. error values Zi-yi for each sample) . 

The error vector z-y and the variance are supplied to the 
Estim ^ block 36. The Estim j8 block 36 also receives the signal 
samples and times y^.t^ from the RF circuitry 22. The Estim ^ 
block can be implemented as a processor executing a program or 
by other suitable means. It operates to generate a Doppler 
compensation factor jS which is a complex number having real and 
imaginary parts* Thus, jS may be represented as: 

Qf+j0t 

where a represents amplitude variations due to fading and 0 is 
the Doppler angle. The Estim jS block 36 also includes a memory 
for holding values required to generate the Doppler correction 
factor 0, including C^, B{0) and E{2-y) . These terms are 
explained later with reference to Figure 3 . 

The Doppler correction factor 0 together with the signal samples 
yi and times ti for each burst are supplied to the transforming 
circuitry 40 to implement Doppler correction according to the 
following equation: 

Yn-compensaced " Yn^l^^^n^ • 

Another implementation of the transforming circuitry 40 is to 
separate the amplitude and phase components into separate" 
correction steps. The phase can be corrected by multiplying the 
samples with e"^*^'^^^^^^^, and the amplitude can be corrected in a 
separate step : 

Jf n-phase-compensated Jrn ^ n 

Yn-compensated = Y-phase- compensated { 1 -real ( /3 ) t^) (equation 5) 

According to the above, imag (/3) and real {/?)=Qr. Thus, the 
received signal is either amplified or alternated during the 
burst depending whether real (/S) is negative or positive. 
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This generates a set of compensated signal samples Ycomp which can 
be fed to the filtering/equalising circuitry 12 to generate 
refined values for the data of the burst. 

As an alternative to the above, in which the training sequence 
X is used to supply known symbols as a reference, bits estimated 
in an equalisation circuit 30 can be used. The equalisation 
circuit 30, for example of Viterbi equaliser, receives the 
filtered, demodulated and equalised signal DATA (r) from the 
filtering and equalisation circuitry 12. The equalisation 
circuit 3 0 operates on a part of the data sequence data of the 
burst (that part having been derived from Estim. block in Figure 
1) to estimate output bits which were sent from the mobile 
station. This output is referred to herein as estim_bits. The 
equalisation circuit 30 operates to make decisions of the bits 
as in known mobile coramunication systems and thus it will not be 
described further herein. These estimated bits can be used in 
place of the training sequence to supply known values for 
calculating the reference vector z. The estim_bits are supplied 
to a reference circuit 34 with the channel impulse response h 
which generates an error vector z'-y, where z '•=estim_bits*h. 

The theory underlying the generation of the Doppler correction 
factor iS will now be outlined. 

The problem of (pure) Doppler frequency estimation can be written 
as : 



Xh=Yr+w, where r is a rotation vector, 
Y is a diagonal sample matrix, and 
w is noise. 



Yo 
0 



Y= 



0 

Yi 



0 

0 

Yn-i 



0 



n-1 
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Xh corresponds to convolution between known symbols {diagonal 
matrix x of training sequence symbols) and impulse response h. 

So, in the Doppler frequency estimation, the Doppler angle 0 is 
to be estimated. This problem is nonlinear and therefore a bit 
tricky to be solved, particularly for implementation in real 
time. A linearisation step can be perfoirmed assuming e^^*=Q«l+j<^to. 
Due to fading the required Doppler frequency correction cannot 
be assumed to be constant in .amplitude. 

If the term j^to has also a real component, the amplitude 
variations can roughly estimated. Thus, the term j^tg is 
replaced by a more general assumption l+^tg, where 0 is complex 
(af4.j<^t) . 

The linear problem can be written as (for simplification 

2=Xh= [Zo^i- . - 

Zi=yi(l+iSti)+Wi, 
which can be reformatted in matrix notation as 

z-y=pj8+w (equation 6) 

where p= [yotoyxti . . .yn-it^-i] ^ and jS is to be estimated. 
A solution to equation 6 can be found by using an LMMSE estimator 

where E(i8) is an expected value for jS, C^^.y), Ct^.yj^^.y) are 
covariance matrixes between the two elements, viz E{j8(z-y)^) and 
E((z-y> (z-y)^ and E(z-y} is the expected error in the received 
burst . 

In this phase it can be assumed that the noise is white gauss ian 
noise having a variance of 5^ and after some manipulation 
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jS=E(^)+p^(z-y-E(2-y) )/(5VC^+p'^) . where 
C.=E( (^-E{j8) ) {/3-E{^) )«) 



(equation 7) 
(equation 8) 



is the complex conjugate transpose matrix of p. 

Equation 7 is used in the Sstim jS block 36 to generate /S. The 
Estim block is illustrated schematically in Figure 3 . In 
practice^ it can be implemented as software executed by a 
processor with a memory foi; holding the required data. An 
Estimate_j8 block 50 operates on the data (Yi^ti) received in each 
time slots to generate a value for jS which is used by the 
transforming circuitry 40, An Acc ^ block 52 accumulates and 
averages ^ values over preceding time slots to generate E(^) . 
The Acc B function can be implemented by an adaptive filter. An 
Acc (2-y) block 54 generates an expected error E(z-y). To do 
this, it can use values for z-y from preceding time slots, or 
E(/3} from the Acc jS block 52 to estimate a frequency error. A 
Gen block 56 produces a value for C^. A number of 

possibilities exist for the parameters of equation 7* 

If E{0) is set to zero, and is also set to zero^ equation 

7 resolves to a normal linear regression. 

As mentioned above, it is intended that ^ will be generated for 
each time slot. The expected value Ei/S) can be generated by 
calculating the average of estimated ]S values from previous time 
slots. Thus, these values can be retained as they are calculated 
for each time slot, and used to produce an average value held in 
the B(/3) block in the Estim jS block 36. 

Another simplification is to assume that the real part of 

is zero. This is a reasonable assumption because ; in typical 

fading channel situations, no long time amplitude is expected. 

Thus, it is adequate to perform long time averages only in 

respect of error components of frequency. Also, the real part 

could be assumed to be zero if only phase variations are to be 

corrected. 
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It is possible to avoid having to calculate the expected error 
E(z-y) by derotating the received samples by e"^*^^^^^'^. For 
this purpose, can be an earlier estimated value of (3 from a 
previous time slot. In this case, equation 7 changes to 

i8=E +p« (z-y) / (5VC^+P^) (equation 9) 

Thus, this requires a preprocessing step (see 50a in Figure 3) 
in which the sample values y^. are derotated by e'^*^™"^*^^^ prior to 
forming the vector p. 

The value for C(j8) can be generated by the GEN block using 
different approaches. It can be a constant value. I can be 
estimated using the information of previous time slots, using for 
example equation 8, It can be made a function of the expected 
value of jS and/or the noise variance 5^. It can be estimated 
using an adaptive filter, or any combination of the preceding. 

The technique is particularly suitable for short delay spread 
channels and pure frequency errors. 
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CLAIMS: 

A -system for generating a Doppler correction factor for 

compensating for Doppler shift in a signal transmitted between 

a mobile station and a base station in a mobile communication 

system, the system comprising: 

a channel impulse response detennnination circuit for 

determining a channel impulse , response for the channel on which 

the signal is received; 

reference circuitry for generating a reference vector from 

a set of known symbols amd the channel impulse responses- 
means for generating an error vector from the reference 

vector and received samples of the signal; 

means for determining noise variance of the incoming signal; 

and 

an estimator for generating a Doppler correction factor 
having a real part representing amplitude correction and an 
imaginary part representing phase correction, the estimator 
utilising the noise variance, the error vector, the samples of 
. the received signal and the sampling times. 

2. A system according to claim 1, wherein the Doppler 
correction factor 13 is generated according to the following 
equation : 

j8=E (iS) -Hp« (z-y-E (z-y) ) / (SVC^+P^^) 

wherein E((3) is an expected value of 

is a complex conjugate transpose matrix of p, where p is the 
vector (yi/ti)*^, 

2-y is the error vector, 

E{2-y) is the expected error in the received samples, 
5^ is the noise variance, and 
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is a constant , 

3. A system according to claim 2, wherein E{I3) and 5^/C^are set 
to zero. 

4. A system according to claim 1, 2 or 3 , which comprises means 
for applying a Doppler correction to the received signal samples 
using the Doppler correction factor jS. 

5. A system according to any preceding claim wherein the signal 
comprises a transmission burst in a time slot in a TDMA mobile 
communication system, and wherein a new Doppler correction factor 
jS is generated for each time slot. 

6. A system according to claim 5, which comprises means for 
estimating Ei^) as an average of generated values for jS from 
previous time slots . 

7. A system according to claim 5, wherein the real part of E{0) 
"is set to be zero* 

8. A system according to any preceding claim, which comprises 
means for derotating the received signal samples by e'^*^"^^^^*^*^^, 
where E{^) is an expected value for the derotated received 
signal samples then being used to generate the Doppler correction 
factor iS in the present time slot according to the following 
equation: 

^=:E{i8)-^p«(z-y)/(aVC^+P^) . 

S. A method for compensating said Doppler shift in a signal 
transmitted between a mobile station and a base station in a 
mobile communication system, the method comprising: 

determining a channel impulse response for the channel on 
which the signal is received ; 

generating a reference vector from a set of known symbols 
in the channel impulse response; 
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generating an error vector from the reference vector and 
received samples of the signals; 

determining noise variance of the incoming signals; 

estimating a Doppler correction factor having a real part 
representing amplitude correction and an imaginary part 
representing phase correction, the estimating step using the 
noise variance, the error vector, the samples of the received 
signals and the sampling times; and 

applying the Doppler correction factor to the received 
signal samples thereby to compensate for Doppler shifts 

10. A method according to claim 9, wherein the signal comprises 
a sequence of transmission bursts in a TDMA mobile communication 
system, and the Doppler correction factor is generated for each 
transmission burst. 

11. A method according to claim 9 or 10, wherein the Doppler 
correction factor 0 is generated according to the following 
equation: 

je=.E O) +p« {z-y-E (2-y) ) / (dVc^+p"^) 

wherein E(i8) is an expected value of /S, 

is a complex conjugate transpose matrix of p, where p is the 
vector (yi#ti)^, 

2-y is the error vector, 

E{2-y) is the expected error in the received samples, 
6^ is the noise variance, and 
Cfl is a constant. 
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